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Abstract 

A lower bound on Am 2 for atmospheric neutrino oscillations can be obtained 
from the global value of R at a given energy. This bound may be more reliable 
than one based on the observed isotropy. Some global values of R imply an 
anisotropy in the neutrino flux that may or may not be manifest in the ob- 
served neutrino interactions. A limit of Am 2 > 0.0077 eV 2 is easily obtained. 
Much higher mass scales are implied by the reported energy independence 
and magnitude of R 

Subject headings: Cosmic Rays — Elementary Particles — Neutrino Os- 
cillations 
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The recent announcement of the discovery of neutrino oscillations by the Super- 
Kamiokande collaboration Q does not, under closer scrutiny, seem to be consistent with 
the observations. The oscillation length they infer corresponds to 1000's of kilometers at 
their energies. It has been pointed out || that since the oscillation length is large compared 
to half the distance scales in the problem and small compared with the other half that a 
serious up-down anisotropy would be expected in the data. Such a large anisotropy is not 
observed. 

It has been suggested that the apparent near isotropy may be an experimental artifact, 
a consequence of the scattering angle between the neutrino direction and the reconstructed 
lepton, compounded by experimental resolution. While such angular resolution effects should 
decrease with energy (which unfortunately couples the neutrino mass scale to resolution 
effects) in this paper we explore what can be learned if the directional information is ignored. 

We employ the variable R commonly defined to be: 

R = (n/e) DA TA 

(/V e W 

and often interpreted as the fraction of muon neutrinos observed relative the expected num- 
ber. A value of R less than one is often taken as evidence for neutrino oscillations and is 
modeled by: 

P v ^ Vt = sin 2 (2#)sin 2 (1.27Am 2 ^) 

sin 2 (26 l ) and Am 2 are the two parameters of the model. L is the distance, in kilometers the 
neutrino has traveled, E is the neutrino energy, in GeV and Am 2 is measured in eV 2 . 

Only projections of R as a function of energy and of the cosine of the zenith angle have 
been published. It is possible that the observed isotropy is a consequence of the large event 
rate at low energies where the reconstructed track direction is less well correlated with the 
initial neutrino direction. This paper points out that isotropy of the neutrino flux within 
the neutrino oscillation hypothesis is a consequence of the small value of R observed and 
the bimodal nature of the path length scales in the problem. Once R gets below about 
0.73 the upper hemisphere (downward going) flux must participate in the oscillations. Due 
to the narrow range of distances for the downward component most of the solid angle will 
participate in oscillations. 

Underground detectors can observe atmospheric neutrinos. For atmospheric neutrinos 
the distance the neutrino has traveled depends on the direction it enters the detector. Since 
the detectors are located near the surface of the Earth path lengths vary from a few kilome- 
ters to about 13,000 kilometers. Figure [I] is a plot of the log of the path length as a function 
of neutrino direction for neutrinos originating in the atmosphere. The path length can be 
approximated by: 

L(cos(# 2 )) = ylR\{cos 2 {6 z )-\) + B% - R x cos(0 2 ) 

with R2 representing the distance from the center of the earth to the upper atmosphere where 
the neutrinos are born and R\ representing the distance from the center of the earth to the 
detector (R2 > R\). Figure Q shows that atmospheric neutrino oscillations are dominated 
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by two distance scales. Neutrinos from above have path lengths of 10's of kilometers. Those 
coming from below have path lengths of the order of 1000's of kilometers. 

Neutrinos that have converted from muon type to tau type or a sterile type would not be 
observed in the detector. This reduction of the flux of muon neutrinos, relative to electron 
neutrinos would be manifest by a value of R < 1. Low values of R have been observed [|3],[4|]. 

One can obtain more information about the nature of the oscillations by measuring R 
in specific energy and/or directional bins. But at very low energies the direction of the 
reconstructed event may not be closely associated with the direction of the initial neutrino 
and so one may be averaging over substantial distances. Distance varies slowly with angle 
near the vertical so this is not necessarily a serious problem. 

The global value of R (for an isotropic atmospheric source) is given by 

1 t 1 o o o \/Ri(cos 2 (6 z ) -l) + Rl-R l cos(6 z ) 
R=l-=-x / sin 2 (2#) sin 2 (1.27Am 2 ^-^ — ^)dcos(^) 

2 7-i E 

Which is a function of Am 2 /E 

Since sin 2 (2#) < 1 one can get a bound on Am 2 /E from the observed value of R at a 
given energy. Figure |2] plots si *iT^fg) as a function of Am 2 /E when the flux is integrated 

over all solid angles. ^jrMg\ is an upper bound on 1 — R. At low values of Am 2 / 'E the 
neutrinos have not had time to oscillate so 1 — R is zero. At very large values of Am 2 /E the 
neutrinos traveling over all paths in the problem have gone through several oscillations so 
that ^rMg\ is |. The interesting structure near the central region of the plot comes about as 
the oscillation length becomes comparable to various distances in the problem. The plateau 
at about 0.25 comes about for that range of oscillation lengths for which the downward path 
is too short for significant changes to have occurred but where the upward direction is large 
compared to the oscillation length. A value of R in the intermediate range of 0.5 < R < 0.73 
implies a flux isotropy since the upper hemisphere must participate in oscillations to get R 
this small. One can still get isotropy with a large R for small values of sin 2 (29). 

The observed value of R [| for all contained events is 0.61±0.03±0.05 (R < 0.68 at 90% 
confidence limit). For this work we need the value of R at a specific energy. From figure 4 
of reference || it can be found that R = 0.63 ± 0.05 at 1.2 GeV. This value gives an upper 
bound of R < 0.69 at 90% confidence level. As seen from figure ^ for this value of R that 
\og 10 (Am 2 / E) is greater than about -11.19 in units of eV or -2.19 in units of eV 2 /GeV. At 
1.2 GeV this gives a bound of Am 2 > lO" 2 - 19 eV 2 /GeV xl.2 GeV > 7.7 x lO" 3 eV 2 . 

The observed values of R are energy independent over a substantial range |||§||. From 
figure H it can be seen that R is expected to be independent of energy over a substantial 
range. In particular for Am 2 / E < 10~ 13 ' 8 eV R will be energy independent at a value 
near one. For Am 2 / E > 10 -9,8 eV R will be energy independent with a minimum value 
near |. The intermediate range plateau 10 -13 < Am 2 /E < 10~ 12 ' eV would also appear 
to be energy independent with suitable energy averaging. In this intermediate range the 
observations would have a noticeable "neutrino twilight" || with the upper hemisphere 
being considerably brighter than the lower one. 

The 5 points in reference || and the 4 multi GeV points from || may be fit to find 
the best values of sin 2 (26 l ) and Am 2 . Physical values of sin 2 (29) must be bounded in the 
range < sin 2 (26 l ) < 1. A formal fit is difficult since the errors in references f| and |J are 
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dominated by systematic effects and are correlated. So it is difficult to get a true confidence 
level associated with the fit. From the global value of R cited above it can be inferred that 
the common systematic error is on the order of 0.047. Removing this common factor and 
doing a \ 2 minimization with the resulting reduced errors yields a best fit to sin 2 (26) = 0.8 
and Am 2 = 0.24 eV 2 . The results of the fit are the points plotted on the right of figure 0. In 
actual fact there is no best fit since all values above about Am 2 > 0.1 eV 2 give a reasonably 
good fit to the data. 

For reference the points plotted near the center of the figure are the Super Kamioka "best 
fit" to the neutrino oscillation hypothesis, sin 2 (26) = 1 and Am 2 = 2.2 x 10~ 3 . The Super 
Kamioka fit could be improved if the points could be lowered to center on the curve. But 
the sin 2 (2$) = 1 bound has been reached and the points can not be lowered. To get good 
agreement for Am 2 = 2.2 x 10~ 3 would require a value of sin 2 (2#) = 1.5 which is unphysical. 

(Note: It is the same data from Super Kamioka plotted at the center and the right of 
figure ^. They appear in different places since the parameters of the fit sin 2 (26) and Am 2 
have been incorporated into the plot axes.) 

It is difficult to quantify the confidence level associated with this fit of Am 2 to R(E) for 
all 9 points. In an 1MB presentation of the atmospheric neutrino anomaly J| the common 
systematic errors are clearly assigned to the expected values and the measured values are 
illustrated with the truly independent statistical errors. 

Our relation above for R(^ L -) assumed that the source was isotropic. That is, it was 
assumed that the neutrinos were uniformly distributed over the path length distribution of 
figure |l|. At modest energies the neutrino flux may be modulated by geomagnetic effects 
that restrict the number of primary cosmic ray particles that can impact the atmosphere 
and create neutrinos. This would modulate the path length distribution. In the U — D 
analysis of reference [IJ (which is substantially the same data as in |5||)[]) it is indicated that 
such geomagnetic effects are less than 2% at sub-GeV energies and less than 1% at higher 
energies. 

To get agreement with the observed bound of R < 0.68 with a Am 2 of the order of 
2 x 10 -3 eV 2 would require a large neutrino flux anisotropy. An expected upward going 
muon neutrino flux approximately 80% greater than the downward flux (1.8 times) is needed. 
While neutrino oscillations would remove much (50%) of the excess upward muon neutrino 
flux the electron neutrino signal would reflect this serious distortion. To achieve a value of 
R = 0.61 would require 3.5 times greater upward neutrino flux. 

In the absence of very large directional modulation in the parent neutrino flux the upper 
hemisphere needs to participate in neutrino oscillations to achieve such a small value of R 
via this simple two component oscillation mechanism. 
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FIG. 1. Approximate log 10 of the neutrino flight path (in kilometers) as a function of the 

cosine of the zenith angle. 
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FIG. 2. g^iT^fl) plotted as a function of the log 10 of Am 2 /E v with Am 2 in eV 2 and E measured 
in eV. sJ^^g) is the maximum global fraction of muon neutrinos that have transformed into a 
noninteracting species when integrated over all atmospheric neutrino path lengths. The points 
plotted near the center of the figure are from references 5 and 6 using the Am 2 from reference 2. 
The points plotted to the right are our fit to Am 2 using the data of references 5 and 6. 
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